Time-domain Implementation of the Optimal Cross-Correlation Statistic for Stochastic Gravitational-Wave Background Searches in Pulsar Timing Data

Supermassive black hole binaries, cosmic strings, relic gravitational waves from inflation, and first order phase transitions in the early universe are expected to contribute to a stochastic background of gravitational waves in the 10^(-9) Hz-10^(-7) Hz frequency band. Pulsar timing arrays (PTAs) exploit the high precision timing of radio pulsars to detect signals at such frequencies. Here we present a time-domain implementation of the optimal cross-correlation statistic for stochastic background searches in PTA data. Due to the irregular sampling typical of PTA data as well as the use of a timing model to predict the times-of-arrival of radio pulses, time-domain methods are better suited for gravitational wave data analysis of such data. We present a derivation of the optimal cross-correlation statistic starting from the likelihood function, a method to produce simulated stochastic background signals, and a rigorous derivation of the scaling laws for the signal-to-noise ratio of the cross-correlation statistic in the two relevant PTA regimes: the weak signal limit where instrumental noise dominates over the gravitational wave signal at all frequencies, and a second regime where the gravitational wave signal dominates at the lowest frequencies.Comment: 12 pages, 3 figures, submitted to Physical Review

PINT: A Modern Software Package for Pulsar Timing

Over the past few decades, the measurement precision of some pulsar timing experiments has advanced from ~10 μs to ~10 ns, revealing many subtle phenomena. Such high precision demands both careful data handling and sophisticated timing models to avoid systematic error. To achieve these goals, we present PINT (PINT Is Not Tempo3), a high-precision Python pulsar timing data analysis package, which is hosted on GitHub and available on the Python Package Index (PyPI) as pint-pulsar. PINT is well tested, validated, object oriented, and modular, enabling interactive data analysis and providing an extensible and flexible development platform for timing applications. It utilizes well-debugged public Python packages (e.g., the NumPy and Astropy libraries) and modern software development schemes (e.g., version control and efficient development with git and GitHub) and a continually expanding test suite for improved reliability, accuracy, and reproducibility. PINT is developed and implemented without referring to, copying, or transcribing the code from other traditional pulsar timing software packages (e.g., Tempo/Tempo2) and therefore provides a robust tool for cross-checking timing analyses and simulating pulse arrival times. In this paper, we describe the design, use, and validation of PINT, and we compare timing results between it and Tempo and Tempo2

Your: Your Unified Reader

The advancement in signal processing and GPU based systems has enabled new transient detectors at various telescopes to perform much more sensitive searches than their predecessors. Typically the data output from the telescopes is in one of the two commonly used formats: psrfits and Sigproc filterbank. Software developed for transient searches often only works with one of these two formats, limiting their general applicability. Therefore, researchers have to write custom scripts to read/write the data in their format of choice before they can begin any data analysis relevant for their research. \textsc{Your} (Your Unified Reader) is a python-based library that unifies the data processing across multiple commonly used formats. \textsc{Your} implements a user-friendly interface to read and write in the data format of choice. It also generates unified metadata corresponding to the input data file for a quick understanding of observation parameters and provides utilities to perform common data analysis operations. \textsc{Your} also provides several state-of-the-art radio frequency interference mitigation (RFI) algorithms, which can now be used during any stage of data processing (reading, writing, etc.) to filter out artificial signals.Comment: 3 pages, Published in JOSS, Github: https://github.com/thepetabyteproject/you

SEXTANT X-Ray Pulsar Navigation Demonstration: Additional On-Orbit Results

The Station Explorer for X-ray Timing and Navigation Technology (SEXTANT) is a technology demonstration enhancement to the Neutron-star Interior Composition Explorer (NICER) mission, a NASA Astrophysics Explorer Mission of Opportunity to the International Space Station, launched in June of 2017. In late 2017, SEXTANT successfully completed a first demonstration of in-space and autonomous X-ray pulsar navigation (XNAV). This form of navigation relies on processing faint signals from millisecond pulsars-rapidly rotating neutron stars that appear to pulsate in the X-ray band-and could potentially provide a GPS-like navigation capability applicable throughout the solar-system and beyond. In this work, we briefly review prior SEXTANT results and then present new results focusing on: making use of the high- flux but rotationally unstable Crab pulsar, and using XNAV to estimate position, velocity, and time in the presence of an imperfect local clock

Multi-Messenger Gravitational Wave Searches with Pulsar Timing Arrays: Application to 3C66B Using the NANOGrav 11-year Data Set

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.Comment: 14 pages, 6 figures. Accepted by Ap

A Distant Fast Radio Burst Associated with Its Host Galaxy by the Very Large Array

We present the discovery and subarcsecond localization of a new fast radio burst (FRB) by the Karl G. Jansky Very Large Array (VLA) and realfast search system. The FRB was discovered on 2019 June 14 with a dispersion measure of 959 pc cm⁻³. This is the highest DM of any localized FRB and its measured burst fluence of 0.6 Jy ms is less than nearly all other FRBs. The source is not detected to repeat in 15 hr of VLA observing and 153 hr of CHIME/FRB observing. We describe a suite of statistical and data quality tests we used to verify the significance of the event and its localization precision. Follow-up optical/infrared photometry with Keck and Gemini associate the FRB with a pair of galaxies with r ∼ 23 mag. The false-alarm rate for radio transients of this significance that are associated with a host galaxy is roughly 3×10⁻⁴ hr⁻¹. The two putative host galaxies have similar photometric redshifts of z_(phot) ∼ 0.6, but different colors and stellar masses. Comparing the host distance to that implied by the dispersion measure suggests a modest (~ 50 pc/cm⁻³) electron column density associated with the FRB environment or host galaxy/galaxies

The NANOGrav Nine-year Data Set:Mass and Geometric Measurements of Binary Millisecond Pulsars

We analyze 24 binary radio pulsars in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) nine-year data set. We make 14 significant measurements of the Shapiro delay, including new detections in four pulsar-binary systems (PSRs J0613−0200, J2017+0603, J2302+4442, and J2317+1439), and derive estimates of the binary-component masses and orbital inclination for these MSP-binary systems. We find a wide range of binary pulsar masses, with values as low as ${m}_{{\rm{p}}}={1.18}_{-0.09}^{+0.10}\,{M}_{\odot }$ for PSR J1918−0642 and as high as ${m}_{{\rm{p}}}={1.928}_{-0.017}^{+0.017}\,{M}_{\odot }$ for PSR J1614−2230 (both 68.3% credibility). We make an improved measurement of the Shapiro timing delay in the PSR J1918−0642 and J2043+1711 systems, measuring the pulsar mass in the latter system to be ${m}_{{\rm{p}}}={1.41}_{-0.18}^{+0.21}\,{M}_{\odot }$ (68.3% credibility) for the first time. We measure secular variations of one or more orbital elements in many systems, and use these measurements to further constrain our estimates of the pulsar and companion masses whenever possible. In particular, we used the observed Shapiro delay and periastron advance due to relativistic gravity in the PSR J1903+0327 system to derive a pulsar mass of ${m}_{{\rm{p}}}={1.65}_{-0.02}^{+0.02}\,{M}_{\odot }$ (68.3% credibility). We discuss the implications that our mass measurements have on the overall neutron-star mass distribution, and on the "mass/orbital-period" correlation due to extended mass transfer